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ABSTRACT 
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We present the results of an Eulerian adaptive mesh refinement (AMR) hydrody- 
namical and N-body simulation in a ACDM cosmology. The simulation incorporates 
common cooling and heating processes for primordial gas. A specific halo finder has 
been designed and applied in order to extract a sample of galaxy clusters directly 
obtained from the simulation without considering any resimulating scheme. We have 
studied the evolutionary history of the cluster halos, and classified them into three 
categories depending on the merger events they have undergone: major mergers, minor 
mergers, and relaxed clusters. The main properties of each one of these classes and 
the differences among them are discussed. The collisions among galaxy clusters are 
produced naturally by the non-linear evolution in the simulated cosmological volume, 
no controlled collisions have been considered. We pay special attention to discuss the 
role of merger events as a source of feedback and reheating, and their effects on the 
existence of cool cores in galaxy clusters, as well as in the scaling relations. 

Key words: hydrodynamics - methods: numerical - galaxy clusters - large-scale 
structure of Universe - Cosmology 



1 INTRODUCTION 

Galaxy clusters are crucial pieces in our understanding of 
■ the Universe. They represent the transition between two dif- 
ferent regimes. The first one involves the evolution of per- 
turbations on cosmological scales where the physics of the 
problem is only driven by gravity on the dark matter com- 
, ponent. The second one, is the scenario of galactic scales, 
1 where the action of gravity is combined with the effects of a 
complex gas dynamics and all sort of astrophysical phenom- 
ena (cooling, star formation, feedback, ...). 

Precisely, due to its nature linking difi'erent scales, the 
galaxy clusters provide us with a powerful tool to constrain 
the cosmological parameters, and at the same time, they 
help us to understand the physical properties of the intra- 
cluster medium (ICM) and its interplay with galaxy for- 
mation and evolution. In this sense, galaxy clusters are ex- 
tremely valuable laboratories where explore the connection 
between cosmological scales and the formation and evolution 
of galaxies. 

Th e simplest m odel explaining the properties of the 
ICM l|Kaiserlll98d ) assumes that gravity is the only re- 
sponsible force determining the evolution of the ICM. In 
this scenario, the gas collapses into the dark matter poten- 
tial wells and then, accretion shocks form moving outwards 
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in the cluster and h eating up the gas until the v irial temper- 
ature of the halo (|Quilis. Ibafiez fc Saej[l998l ). As gravity 
acts on all scales, this model has been known as self-similar. 

This model gives us precise predictions for the prop- 
erties of galaxy clust ers: X-ray lumi nosity, temperature, 
entropy, and mass (iBryan fc Norman 1998). However, 
the scaling relations produced by the self-similar model 
do not match the observational results completely. More 
specifically, the slope of the X-ray lumi nosity-temperature 
relation is steeper tha n predicted (iMa rkcvitch' Il998l : 
lArnaud fc EvraridI Il999l : lOsmond fc Pomnan .20041), the 
measured gas entropy in poor clusters and groups is 
higher than expected (jPonman. Sanderson fc Finoguenovl 
1^03), and it has been observed a decreas ing trend of 
the gas mass fraction in poorer systems (iBaloeh et al.l 
I2OOII : iLin. Mohr fc Stanford! l2003l : ISanderson et al.ll2003l ). 
In addition, there is an important scatter in the relations 
I Fabian et al.|[l993 ) partly, but not totally, connected with 
the effect of the different environments where clusters live. 

Those discrepancies between the self-similar model and 
the observations have motivated the idea that some impor- 
tant physics, basically related with the baryonic component, 
is missing in the model. 

Some of these non-gravitational processes have been 
included in simulations trying to solve the similarity 
breaking: preheating ([Navarro. Frenk fc White! Il995l : 
iBialek. Evrard fc Moh^ I2OOII: iBorgani et al.! 12002!). and 
radiative cooling ( Pearce et al.! 2000l : Muanwong et al.] 
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Merger events can also be an important source of feed- 
back in galaxy clusters. They can produce shocks and com- 
pression waves in the halos which eventually can release part 
of the energy associated with the collision as thermal en- 
ergy in the final system (|McCarthv et aL I l2007l ). It is likely 
that turbulence and mixing could play an important role in 
how this energy is mixed and released in the ICM of the 
final halo after the merger. As it is well known, some of 
the results of these simulations could depend on the ability 
of different numerical techniques to describe shock waves, 
strong gradients, turbulence, and mixing, which can be very 
different. Although is still a matter of debate, it has been 
shown, at least for some idealised tests, that the compari- 
son between grid codes and SPH codes - when numerical 
resoluti on is similar - can give substantial dif ferences in the 
results (|Frenk et al.lll999l : lAgertz et al.ll2007l '). It seems rea- 
sonable to think that these inherent numerical differences 
could translate into relevant differences when they are ap- 
plied to more complex and realistic scenarios like galaxy 
clusters. This situation makes interesting, and complemen- 
tary, to pursue the number of studies using the different 
numerical strategies available. 

In this paper, we want to investigate the role of the 
galax;y cluster mergers as a source of feedback and re- 
heating in a complete general cosmological framework. The 
galax;y clusters form and evolve due to the non-linear evo- 
lution of primordial perturbations, and therefore, no spe- 
cial symmetry or idealised clusters are considered. In this 
scenario, the merger events naturally take place according 
to the hierarchical evolution. Previous works have exten- 
sively studied t he mergers of galaxy clus t ers using con trolled 
collis i ons (e.g. 'Ricker fc Sarazin (2001); Poole et al.1 l|2006l . 
120071 ) : lMcCMthv et al. (2007.): , Poole ct al. (20081 )). The ap- 
proach adopted in the present paper could be considered as 
complementary to the studies using controlled mergers. It 
is clear that our approach has some important weaknesses, 
when it is compared with controlled mergers, like the worst 
resolution or the impossibility to control the different param- 
eters involved in the problem. However, it gives a description 
of the problem in a cosmological context, without symme- 
tries, including the presence of substructures and taking into 
account the effects of the different environments. 

In order to study the role of mergers, fulfilling all the 
previous requirements, we have carried out a simulation of a 
moderate size box of side length 100 Mpc. We have per- 
formed it with an Eulerian AMR cosmological code includ- 
ing the usual processes of cooling and heating for primordial 
gas, and a phenomenological star formation treatment. We 
have identified and followed the evolution of the different 
galaxy cluster halos. Once the evolutionary history of the 
halos is known, we have classified them into three broad 
categories depending on the features of the merger events in 
which they have been involved. We will discuss their effects 
on cluster properties. These mergers are the ones naturally 
happening in the building up of the galaxy clusters. 



The paper is organized as follows. In Section 2, we 
present the details of the simulation and describe the halo 
finder used to identify the galaxy cluster halos. In Section 3, 
we analyse the results of the simulation describing the main 
properties of the galaxy cluster sample, and the effects of 
mergers. Finally, in Section 4, we summarise and discuss 
our results. 



2 THE SIMULATION 
2.1 Simulation details 

The simulation described in this paper h as been perf ormed 
with the cosmological code MASCLET (|Quilisll200j ). This 
code couples an Eulerian approach based on high-resolution 
shock capturing techniques for describing the gaseous com- 
ponent, with a multigrid particle mesh N-body scheme for 
evolving the coUisionless component (dark matter). Gas and 
dark matter are coupled by the gravity solver. Both schemes 
benefit of using an adaptive mesh refinement (AMR) strat- 
egy, which permits to gain spatial and temporal resolution. 

The numerical simulation has been performed assum- 
ing a spatially fiat KCDM cosmology, with the following 
cosmological parameters: matter density parameter, f2m = 
0.25; cosmological constant, Oa = A/3H^ — 0.75; baryon 
density parameter, Qb = 0.045; reduced Hubble constant, 
h = Ho/lOOkms-^ Mpc-^ = 0.73; power spectrum index. 
Us — 1; and power spectrum normalisation, erg = 0.8. 

The initial conditions w ere set up at 2: = 50, u sing a 
CDM transfer function from lEisenstein fc HrJ (|l998l ). for a 
cube of comoving side length 100 Mpc. The computa- 
tional domain was discretized with 512^ cubical cells. 

A first level of refinement (level 1) for the AMR scheme 
was set up from the initial conditions by selecting regions 
satisfying certain refining criteria, when evolved - until 
present time - using the Zeldovich approximation. The vol- 
umes selected as refinable were covered by grids (patches) 
with numerical cells selected from the initial conditions. The 
regions of the box not eligible to be refined were degraded in 
resolution by averaging the quantities obtained on the initial 
grid. This procedure creates the coarse grid (level 0) for the 
AMR scheme. These coarse cells have a volume eight times 
larger than the first level ones. In the same manner, the dark 
matter component within the refined regions was sampled 
with dark matter particles eight times lighter than those 
used in regions covered only by the coarse grid. During the 
evolution, regions on the different grids are refined based on 
the local baryonic and dark matter densities. Any cell with a 
baryon mass larger than 5.6 x 10* Mq or a dark matter mass 
larger than 2.5 x lO^Af0 was labelled as refinable. The ratio 
between the cell sizes for a given level {I + 1) and its parent 
level (l) is, in our AMR implementation, Axi+i/ Axi — 1/2. 
This is a compromise value between the gain in resolution 
and possible numerical instabilities. This method produces 
patches with a boxy geometry and cubic cells at any level. 

The simulation presented in this paper has used a max- 
imum of seven levels (/ = 7) of refinement, which gives a 
peak spatial resolution of 3 kpc. For the dark matter we 
consider two particles species, which correspond to the par- 
ticles on the coarse grid and the particles within the first 
level of refinement at the initial conditions. The best mass 
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Table 1. Cluster sample. Main features of selected clusters at 2 = 0. Column 2, virial radius in 
units of Mpc] column 3, total mass within the virial radius in units of lO^* Mq; column 4, 
mass-weighted temperature within the virial radius in keV; column 5, gaseous mass within the 
virial radius in units of lOi* A/q; column 6, average entropy within the virial radius in units of 
keVcm?; column 7, type classification: relaxed or with no important mergers (R), minor mergers 
(MI), and major mergers (MA). 



cluster 


(/i~i Mpc) 


l^iU ^^-'-01 


T 

(keV) 




S 


type 


CLOl 


2.32 


18.61 


7.02 


3.45 


2601.06 


MA 


CL02 


2.22 


15.83 


5.99 


4.89 


2130.17 


MA 


CL03 


1.58 


5.68 


3.24 


1.08 


919.34 


MA 


CL04 


1.48 


4.70 


3.80 


0.79 


998.65 


MA 


CL05 


1.39 


3.93 


2.29 


0.70 


1063.19 


MA 


CL06 


1.07 


1.86 


1.26 


0.37 


510.01 


MA 


CL07 


1.01 


1.53 


1.01 


0.25 


846.33 


MA 


CL08 


0.93 


1.14 


0.87 


0.21 


422.39 


MA 


CL09 


1.51 


5.19 


3.09 


1.03 


1166.28 


MI 


CLIO 


1.51 


5.11 


2.18 


0.92 


1126.22 


MI 


CLll 


1.36 


3.76 


2.55 


0.49 


862.07 


MI 


CL12 


1.45 


4.40 


2.52 


0.97 


980.39 


R 


CL13 


1.10 


1.99 


1.40 


0.41 


409.21 


R 


CL14 


0.99 


1.39 


1.18 


0.29 


362.36 


R 


CL15 


0.89 


1.02 


0.97 


0.22 


554.70 


R 


CL16 


0.89 


1.01 


0.90 


0.22 


386.56 


R 



resolution is 5.75 x 10* h ^ Mq, equivalent to distribute 512^ 
particles in the whole box. 



Our simulation includes cooling and heating processes 
which ta ke into account Compton and free-free cooling, UV 
heating (jHaardt fc MadaulirOQ^ ). and atomic and molec- 
ular cooling for a primordial gas. In order to compute 
the abundances of each species, we assume that the gas 
is optically thin and i n ionization equilibrium, but not 
in thermal equilibri um (|Katz. Weinberg fc Hernauist|[l99^ : 



pTh euns et ahl Il998f). The tabulated c ooling rates were 
taken from " iSutherland fc Dopital (jl993l ) assuming a con- 
stant metallicity 0.3 relative to solar. The cooling curve was 
truncated below temperatures of 10* K. The cooling and 
heatirig were included in the energy equation (see Eq.(3) in 
iQuilisI l|2004 )) as extra source terms. 



The star formation has also been modelled with the 
phenomenological app roach c o mmon l y used in cosmolog ical 
simula tions (i.e. Yepe s et al. I l|l997l ): ISpringel fc Hcrng uistl 
l|2003l )). A more detailed description of how star formation 
has been introduced in MASCLET code is presented in Ap- 
pendix A. Despite the use of an AMR code to perform the 
simulation described in the present paper, we have still had 
numerical limitations, namely, the number of patches placed 
at the highest level of refinement. Although this limitation 
has been not crucial for the description of clusters, it has 
translated into a poor star formation efficiency as the anal- 
ysed run only allowed star formation at the highest level 
of refinement. This apparent drawback of our simulation is 
not dramatic for the purpose of the present study where we 
focus in the effect of mergers on the ICM properties. Thus, 
the stellar feedback has turned out to be very low, and con- 
sequently, it does not alter the pure effect of the mergers. 



2.2 Cluster identification 

A crucial issue in the analysis of our simulation has to do 
with the cluster identification. In order to do so, we have 
developed a halo finding method specially suited for the fea- 
tures of the cosmological code MASCLET. 

The halo finder developed for M ASCLET code foUows 
a sim ilar idea to MHF halo finder (|Gill. Knebe fc GibsonI 
12004'). The main aim is to identify gravitationally bound 
objects in a N-body simulation. 

We have used an identification technique based on the 
original idea of the spherical over-density (SO) method. The 
basic concept of this technique is to identify spherical regions 
with an average density above a certain threshold, which can 
be fixed according to the spherical top-hat collapse model. 
Therefore, we can define the viral mass of a halo, Al^ir, as 
the mass enclosed in a spherical region of radius, r^ir, having 
an average density Ac times the critical density pc{z) = 
3H{zf/8nG: 



4 3 

Mvir{< rvir) — -TTT^irAcPc- 



(1) 



The over-density Ac depends on the adopted cosmo- 
logical mo del, and can be approx imated by the following 



logical mo del, and can be approx in 
expression (|Brvan fc Norman|[l998l '): 



Ac = 18tv^ + 82x - 39x^ 



(2) 



where x = n{z) - 1 and n(z) = [0,m{l + z)^]/[r2m(l -I- 
z)^ -I- ^a]. Typical values of Ac are between 100 and 500, 
in particular, for the cosmological parameters considered in 
our simulation, Ac — 373. 

The SO method implies oversimplifications which could 
lead to somehow artificial results which deserve a careful 
treatment. The enforced spherical symmetry for halos that, 
at least at some mome nts of their liv es, can have quite ir- 
regular structures (e.g., IWhitd (|2002l )). or the difficulties to 
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MAJOR MERGER (CL03) 



MINOR MEROER (OLIO) 



RELAXED HALO (0L1 2) 



[0-25[ % _ 
[25-50[ % _ 
[50-75[ % .. 
[75-100] %_ 



(5 ® o 




Figure 1. Merger tree for one cluster of each category: left panel shows a major merger cluster, central panel represents a minor merger 
cluster, and right panel stands for a relaxed cluster. Cluster halos are represented by circles whose sizes are normalised to the final mass 
at z = 0. Lines connecting halos at different times indicate the amount of mass transferred from the progenitors to their descendants. 



discriminate among close density peaks, are some of these 
well known problems. 

The practical implementation of our halo finder has 
several steps designed to improve the performance of the 
SO method and to get rid of the possible drawbacks. In a 
first step, the algorithm finds halos by centring at density 
peaks and growing spheres until the average over-density 
falls below Ac, or there is a rising up of the slope of the 
density profile. The second step takes care of possible over- 
laps among the preliminary halos found in the first step. 
In our method, halos which overlap more than 30% and 
less than 80% in mass are joined as a single one. Halos 
sharing more than 80% of their mass are considered as a 
misidentification and one of them drops out of the list. The 
third step checks that all particles contained in a halo are 
bound. In order to determine whether a particle is bound 
or not, we estimate the escape velocity at the position of 
the particle (jKlvpin. Gottlober fc Kravtsovll 19991 ). If the ve- 



locity of a particle is larger than the escape velocity, the 
particle is assumed to be unbound. A final step verifies 
that the density profile is con sistent with a NFW profile 
iNavarro. Frenk fc Whit3l 19971 '). 

One of the main advantages of our method is that the 
structure of nested grids created by the AMR scheme al- 
ready follows the density peaks, and therefore, densities are 
already calculated by the cosmological code. Other impor- 
tant point, inherent to the structure of the AMR scheme 
used, is that no linking length is needed. The process of 
halo finding can be performed, independently, at each level 
of refinement of the simulation. Then, in a natural way, our 
halo finder can trace halos-in-halos and obtain a hierarchy 
of nested halos. 

The different progenitors are identified by following all 
particles belonging to a given halo backwards in time. This 
procedure can be repeated until the first progenitor of a 
certain halo is found. This method allows us, not only to 
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know all the progenitors of each halo, but the amount of 
mass received from each one of its ancestors. 



3 RESULTS 

In our simulation, we have identified more than three hun- 
dred galaxy clusters and groups spanning a range of masses 
from 1.0 X 10^^ Mq to 2.0 x 10^^ Mq. We refer to this set of 
clusters as the complete sample. We have constructed their 
evolutionary histories and, based on their merging histories, 
we have classified them into three categories according to 
the mass ratio of the halos involved in the collision. 

It is also convenient to adopt a timescale limit since 
mergers occurring at a very early epoch would not have any 
important consequence on the present properties of the clus- 
ters. In this sense, and only for the purpose of delimiting the 
merger events happening recently, we have defined the for- 
mation redshift of a cluster, Zfor, as the redshi ft at which the 
clust er mass is half of its present virial mass (|Lacev Sz Cold 
1 19931 ). Thus, we consider for each cluster only those mergers 
that have relevant effects on its recent past. 

Therefore, taking into account the formation redshift of 
the clusters, Zfar, and the masses of the most (less) massive 
halo, m2 (mi), involved in the merger, we have classified the 
clusters into three categories: 

• Major mergers. Those systems where the mass ratio 
7712 : mi is smaller than 3 : 1. Therefore, a major merger 
involves clusters with similar masses. 

• Minor mergers. Those systems where the mass ratio is 
3 : 1 < m2 : mi < 10 : 1. 

• Relaxed halos. Those systems which have suffered merg- 
ers with very small halos, 10 : 1 < m2 : mi, or smooth 
accretion. 

Out of the complete sample, we have picked up a sub- 
sample which contains the sixteen most massive galaxy clus- 
ters in the computational box. They constitute what it 
would be referred to as the reduced sample, and their main 
features are summarised in Table [1] Depending on the par- 
ticular analysis, that we will be interested in the following 
sections of this paper, we will use the complete or the re- 
duced sample, respectively. Concerning their merging clas- 
sification, in the reduced sample we have found five relaxed 
clusters (R), three have been categorised as clusters with 
minor mergers (MI), and eight have been classified as major 
merger systems (MA). 

Despite we have used an AMR code to carry out the 
simulation described in this paper, due to numerical limita- 
tions, a biased sample of clusters - with a tendency to better 
describe the most massive ones - has been produced. These 
artificial results could be overcome, in future applications, 
by performing resimulations of the selected clusters in the 
sample, although this could prevent us from following the 
mergers in a cosmological context. In any case, we consider 
that the effects of mergers would be more important in those 
systems with higher masses - well described in the present 
simulation - and therefore, we believe that this bias has only 
minor consequences. 

In order to analyse the results of our simulation, we 
will study several thermodynamical properties which can be 
directly connected with observational data, and which have 



been widely studied by all sort of different simulations. In 
addition to the common plots of density, we will also study 
the behaviour of the ICM temperature. X-ray luminosity, 
entropy, and the internal and kinetic energies. 
The ICM temperature will be defined as: 



T = 



(3) 



where Ti and Wi are the temperature and the weight given 
to each cell. In most of the applications in the present pa- 
per, the weight will be the cell mass, Wi — rrii, and therefore, 
this will be a mass-weighted temperature. In some particu- 
lar cases, and for the sake of comparison with observational 
data, we will also use the so cal led spectroscopic-like tem- 
perature (|Mazzotta et al.ll2004l ). Tgi, where the weight is 
Wi = mipiT~^^'^ with pi the density at the cell i. 

A crucial observable quantity, directly related with the 
temperature and the density of the gas, is the bolometric 
X-ray luminosity. In simulations, this quantity can be com- 
puted by adding up all the contributions from each elemental 
volume of gas: 



/ n,niA{T,Z)dV 
Jv 



(4) 



where rie and rii are the electron and ion density, re- 
spectively, and A is the normalised cooling function de- 
pending on the temperatur e (T) and metallicity (Z) from 
[Sutherland fc Dopital (|l993l ) (see Sec. 2.1). 

The next thermodynamical quantity we will pay a spe- 
cial attention is the entropy, which is a n extremely useful 
quantity providing a lot of information (| Void [20051 ) about 
the evolutionary state of the clusters, since it records the 
thermodynamical history of the ICM produced by the grav- 
itational and non-gravitational processes. We will adopt the 
following common definition for the entropy: 

5' = fcflTr^J^/^ (5) 

where rie is the electron number density and fcs is the Boltz- 
mann constant. 

Other thermodynamical quantities useful to quantify 
the effects of mergers and shocks as a source of feedback 
are the total internal energy, Eu, and the total kinetic en- 
ergy, Ek, which are given by the following expressions: 



/ pedV 



Ek 



pv dV 



(6) 



(7) 



where p, e and v are the gas density, specific internal energy 
and the peculiar velocity of the gas fiuid element, respec- 
tively. 

For the sake of comparison with observational data, it 
is also useful to define some characteristic quantities widely 
used in the literature. Similar to the definition of the virial 
quantities (see Sec. 2.2), we introduce a characteristic radius, 
r^, such that the mean density enclosed within this radius 
is A times pc, and therefore, the mass is: 

4 



Ma{< 



-Tvr^Apc. 



(8) 



Consequently, we follow the common definition for the 
temperature. 






Figure 2. Average radial profiles for all clusters belonging to each class: major mergers (left column), minor mergers (central column), 
and relaxed (right column). From top to bottom, the first and second rows display gas (p/ P200) ^^"^ dark matter {pdml P^oo) densities, 
respectively, the third row shows mass- weighted temperature {T/T^qq), and the fourth row represents entropy {S/S2qq)- All profiles have 
been scaled by the mean value of the plotted quantities within the scale radius, t^qq ■ Continuous lines show the average for all the 
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T — 



2r, 



kE 



(9) 



where /i and nip are the mean atomic weight and the proton 
mass. In the same manner, we define the entropy, 



Iinip 

fbpA 



2/3 



(10) 



with ft, = Qt/^m- 

In the following sections, and in order to compare with 
previous works, we will consider A = 200 or A = 500, de- 
pending on the particular case we compare with. All the 
quantities we have just described, are going to be used to 
analyse the results of the simulation. 



3.1 Merger history of selected clusters 

We have constructed the merger tree of the selected galaxy 
clusters tracking all the dark matter particles that belong 
to a given cluster backwards in time. Figure [1] displays 
the merger trees of three halos, which could be considered 
as prototypical ones of each category (i.e., relaxed, minor 
merger, and major merger). The clusters have been selected 
such that they have very similar masses and sizes ai z — 0. 
The merger trees start ai z — and plot all the parent halos 
of the final ones in previous time steps over several output 
times of the simulation. 

The total mass of each halo is represented by a circle, 
whose size is normalised to the mass of the final halo at 
z = 0. The lines connecting circles between diflFerent times 
inform us about the progenitors of a given halo at a given 
time. In addition, the type of line tells us the amount of 
mass transferred from the progenitor to the halo at the con- 
sidered time. Thus, a halo at a certain time connected with a 
progenitor halo at earlier time by a dash-dotted line, means 
that up to 25% of its mass is due to the contribution of that 
progenitor. The same idea follows for other line types. The 
aim of this kind of plot is to show the merger history and 
the different interconnections over time. The horizontal axis 
is designed to separate halos for plotting purpose only, and 
it has no direct implication on the position of halos in real 
space. Vertical axes show the redshift. 

In Figure [T] the diflFerent merger events can be easily 
identified. Whereas the relaxed cluster (right panel) has a 
quiet evolution, the middle panel shows a cluster suffering 
three mergers between z = 0.79 and z = 0.65. By comparing 
the masses of the different halos involved in these processes, 
all the events are classified as minor mergers. In the left 
panel, a cluster undergoing several mergers at different times 
is presented. Some of the mergers are minor ones, but there 
are major collisions ai z = 0.99,2 = 0.79 and z — 0.72. 
In order to quantify the effect of mergers, we will correlate 
all these phases of activity in the clusters evolution with 
changes and effects on the different physical quantities. 

It is important to notice that some of the merger trees 
also show halos that break apart, that is, lose mass and re- 
duce their sizes. This process operates at two well separated 
regimes with different causes. The first group is formed by 
very small size halos. These halos are not really gravitation- 
ally bound and they can be easily disrupted by interactions 
with environment or with other halos. For the halos with 



1.4 
1.2 
1.0 
0.8 
0.6 
0.4 

0.2 

0.0 
0.001 



0.010 



0.100 



1.000 



10,000 



Figure 3. Mean spectroscopic-like temperature profile for all the 
clusters in Table 1. Each cluster profile is normalised by its mean 
spectroscopic-like temperature (Tx ) within t^qq ■ The normalised 
profiles are averaged in order to obtain the plotted mean profile. 
Dots represent the different bins equispaced in logarithmic scale 
with width 0.1 dex, and error bars are 1 a s.d. The shaded re- 
gion encloses the mean and 1 a s.d. tempera ture profiles from a 
representative sample of nearby clusters by iPratt et al.l ||2007|) . 
The small black dots corres pond to the temperatu re profiles of 
the clusters in the sample of IVikhlinin et al. 
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Figure 4. Mean entropy profile for all clusters in T able 1 com- 
pared with the fitting in IVoit. Kay fc BrvanI . The radial 
profiles for each cluster are obtained using equispaced logarith- 
mic bins. All of them are normalised to their respective values 
of S^oo , and then, the mean profile is obtained. Continuous line 
represents the mean pro file and error bars show 1 a s.d. Dashed 
line shows the fitting by IVoit. Kay fc BrvanI ||2005|) . 
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Figure 5. Gas entropy at O.lrjgo as a function of temper- 
ature, T200- Simulated clusters are represented by filled cir- 
cles, wlier eas obse rvational data w i tli er ror bars stand for the 
iPonmanTS andcrso n &: Finoguenovl ||2003| ) sample. In order to 
compare simulations and observations, Tjqq = Tx lias been as- 
sumed. 



Icirger masses, those mass losses are small, and they are as- 
sociated with tidal interactions. 



3.2 Average radial profiles 

In order to analyse the main properties of the simulated 
galaxy clusters, we compute radial profiles for several quan- 
tities. These profiles are centred at the centre of mass of 
each halo and run outwards from the centre to a distance 
slightly larger than the radius, r2oo- The bins are equispaced 
in logarithmic scale with widths 0.1 dex. In all the plots dis- 
playing radial profiles, the radial coordinate is normalised 
to the r2oo at this time. 

In Figure (2] we plot averaged radial profiles for sev- 
eral quantities for the clusters in the reduced sample (see 
Table 1). All the profiles are scaled by the plotted quan- 
tities at r2oo defined according to Eqs.(|81|- (|10|) . The mean 
profiles (continuous lines) are computed by averaging all 
the profiles of the clusters of each class. The right column 
stands for the relaxed clusters, the central column repre- 
sents the minor merger clusters, and the left column dis- 
plays the major merger clusters. The plotted quantities are 
gas (p/p20o) and dark matter {pdm/ p2m) densities, temper- 
ature {T/T200), and entropy {S/S2oo)- The continuous lines 
stand for results at 2 = and shadowed regions mark one a 
deviation. Let us stress that, in Figure[2]and in the following 
ones - unless explicitly stated -, we consider mean profiles 
rather than median profiles. 

A detailed analysis of Fig. [2] shows the main features of 
the three categories in which we have classified the differ- 
ent clusters. The comparison of gas and dark-matter density 
profiles does not show notable difi'erences. Whereas for the 
gas density, the relaxed clusters exhibit a slightly higher 



density at the centre compared with the minor and major 
merger clusters, the behaviour for the dark matter is the 
opposite, having the major merger clusters a higher density. 
In any case, the profiles are consistent with the expected 
characteristics of density profiles for galaxy clusters. 

Concerning the temperature profiles, there are no dra- 
matic differences either. All clusters, in the reduced sam- 
ple, show a central core with an almost constant temper- 
ature and a declining profile outwar ds. This result is com- 
patibl e with observational data (e.g., iDe Grandi fc Molendil 
1200^ ), and with the idea of a quite univer sal temperature 
profile for the galaxy clusters l|Loken et al.l 2002). The ma- 
jor and minor merger cluster profiles have very similar cen- 
tral temperatures, although the isothermal core is larger for 
the major merger clusters. The relaxed clusters have a bigger 
isothermal core with a slightly lower value of the tempera- 
ture compared with the major and minor merger clusters. 

The temperature profiles of the most massive clusters in 
our reduced sample do not exhibit a drop in the temperature 
in their central regions. Apparently, thi s could seem to di ffer 
from the results of the simulations by 'Kav et al.' (2007*) or 
th e observational dat a presented by [Vikhlinin et al. (200^ 
or iPratt et al.l (|2007l ). These last observational results show 
clusters with temperature profiles with drops in their central 
te mperatures. This effect is more outstanding in the case 
of IVikhlinin et"aLl (|2005 l). The nature of this discrepancy 
amid both observational results could be related with the use 
of different instrumentation in order to obtain the data of 
both samples. It must be noticed that wher eas |Pra tt_ et al 



(2007) used XMM-Newton, the results of V ikhlinin et al 
(2005) were obtained using CHANDRA - with higher an- 
gular re solution. Th i s cou ld explain that some cool cores 
in the IPratt et al.l |20o3) sample were not properly re- 
solved. In o rder to compare with the results in iPratt et al.l 
(2007) and IVikhlinin et"aLl (2005), we have calculated the 
spectroscopic-like radial temperature profile for each clus- 
ter of our reduced sample (clusters listed in Table 1) and 
normalised them to their respective mean Tsi within r^aa, 
that we denote as Tx ■ In Figure |3l we compare the mean 
of all these radial temperature profiles, represented by dots 
with error bars (o ne a standard devi ation), with the obser- 
vational results in iPratt et al.] ||2007|) marked as the shaded 
region, and the results in IVikhlinin et af] l|2005l ) represented 
by small dots0. Our results are consistent with these obser- 
vational data i n an average s ense, being slightly more similar 
to the data of iPratt et all (|2007l '). 

We can understand our result, if we keep in mind that 
the clusters we are considering are the most massive ones in 
our sample. As we will discuss in more detail in Sec. 3.4, we 
have found that there is a strong anticorrelation between the 
drop of temperature in the central region and the mass of 
the cluster. Therefore, the larger the mass of the cluster the 
smaller the number of clusters with central gradients of tem- 
perature. However, if all the clusters in the complete sample 
are considered, then a relevant fraction of the population 



It must be mentioned that in IVikhlinin et al.lll2005^ ■ the tem- 
perature profiles are plotted against radial coordinate normalized 
to r-^gg. We have ignored this small correction without relevant 
eff'ects for the purpose of the actual comparison. 
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Figure 6. Radial profiles of entropy (S) and temperature (T) at different redshifts. Left column shows the result for a cluster representing 
the major merger clusters, central column for a minor merger cluster and right column for a relaxed cluster. 



16%) shows temperature profiles with central gradients 
(see Sec. 3.4 for more details). 

More interesting is the analysis of the entropy profiles. 
In all cases, the clusters have entropy cores and profiles 
outside 0.2 r2oo compat ible with a power law S{r) oc r" 
jTozzi fc NormanlbOOlI'l. In previous work carried out by 
IVoit. Kav fc BrvanI (|20o"5l ). the authors performed several 
non-radiative SPH and AMR simulations, and studied the 
main features of the entropy profiles of the galaxy clusters 
in their numerical samples. Besides several differences in 
the inner cores, all clusters in their sample, regardless of 
the numerical technique used, showed very similar entropy 
profiles outside a region around 0.2r2oo- In particular, for 
the AMR simulation, they found that the entropy profile 
in the outer regions can be better fitted by the power law 
S{r) — 1.43 S2oo(r/r2oo)^'^- In Figure 21 we compare the 
mean entr opy profile of all the cluste rs in Table 1 with the 
fitting by IVoit. Kav fc BrvanI l|2005h . Continuous line rep- 
resents our mea n entropy profile wit h la er ror bars. Dashed 
line stands for IVoit. Kav fc BrvanI (|2005| ) fitting. Our re- 
sults seem to be compatible with the fitting in the outer 
part of the profiles, whereas in the inner region, where cool- 
ing and feedback processes could be relevant, difi'erences are 
expected. 

In order to compare our results with observational 
data, we have looked at the values of the entropy at 



0.1 r,rin, 5n 1 , and compared them with previous data 
by iPonman. S anderson fc Finoguenovl l|2003l ). In Figure |S] 
we plot the observational data (with error bars) by 
IPonman. Sanderson fc Finoguenovl l|2003l ) together with the 
values for the clusters in our reduced sample with temper- 
atures higher that IkeV (filled circles). The points rep- 
resenting simulated clusters match well with the observa- 
tional data apart from three clusters that are marginally 
compatible. These three objects turn out to have some pe- 
culiarities as they have suffered quite recent merger events. 
The values of the entropy at the very center of the simu- 
lated clusters are a lso similar to recent Chandra observa- 
tions (|Morandi fc Ettori 2007 '). Therefore, our results seem 
to be reasonably compatible with the observations taking 
into account all the simplifications and limitations of our 
approach. 

Coming back to the comparison of the results according 
to the merger history of the clusters, the generic shape of the 
entropy profiles does not depend systematically on the mass 
or tem perature of the clusters, in agreement wi th observa- 
tions (|Ponman. Sanderson fc Finoguenovl l2003l ). The sizes 
of the cores are similar in the relaxed and minor merger clus- 
ters and slightly larger in the major merger ones. As it would 
be naively expected, the entropy floor in the relaxed clusters 
is lower than for the minor merger clusters, and this one is 
also lower than for the major merger clusters. Although the 
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differences seem not to be dramatic, they are clearly visi- 
ble in the mean profiles. These differences in the value of 
the entropy in the core, would be a clear consequence of the 
different evolutionary histories of each cluster. 

3.3 Merger effects 

Merger events can produce shocks and compression waves in 
halos. Their effects could be an efficient way to transfer part 
of the gravitational energy associated to the collisions to the 
ICM of the final halo after the mergers. In this picture, the 
role of turbulence and mixing phenomena is crucial as a way 
to redistribute this energy into the ICM. 

The study of these scenarios requires a numerical 
scheme able to tackle with an accurate description of shock 
waves, strong gradients as well as to describe the turbulence 
associated to those violent events. As it has been discussed 
in Sec. 1, the ability of different numerical techniques to 
describe these phenomena is still a matter of debate. 

We focus in this Section on the effects of different 
merger events in the thermodynamical properties of the 
ICM. In order to do so, we select the same three clusters 
than in Sec. 3.1. Each one of them represents one of the 
three groups of clusters, and as it was already mentioned, 
they have been chosen in such a way that they have similar 
masses and sizes at z = 0. 

In order to discuss the effects of different merger events 
in the thermodynamical properties of the ICM, we show in 
Figure |6] the radial profiles of entropy {S) and temperature 
(r) at several redshifts for the selected clusters. This figure 
can be correlated with Fig. [1] to detect the effect produced 
by the merger events. Lines representing high redshifts must 
be taken carefully. They correspond to early stages of the 
clusters formation when these structures are far from be- 
ing relaxed and, therefore, the radial profiles are not really 
meaningful. 

The relaxed cluster shows a higher entropy and temper- 
ature at high redshifts, with a tendency to reach a relaxed 
state around z ~ 0.5 with small changes. The tendency for 
the minor merger cluster is similar for the lines displaying 
z ~ 2 and z ^ 1, that is, a reduction of the value of the 
entropy core. However, associated with the minor merger 
events, there is a significant increase in the value of the en- 
tropy core which, eventually, ends up in a reduction of the 
temperature and entropy at z ~ with respect to the values 
at z ~ 0.5. In the case of the major merger cluster, between 
2 ~ 1 and z ~ 0.6 when the major mergers take place, there 
is an increase in the entropy and a reheating as the temper- 
ature also increases. Later on, the cluster cools to a lower 
temperature but part of the energy of the merger has been 
released in the cluster which has a higher entropy. It must 
be noticed that the values for the entropy are considerably 
larger than for the other two clusters previously discussed. 

So as to assess more clearly the effects of mergers in 
the thermodynamical properties of ICM, we show in Fig- 
ure [7] the time evolution of several quantities: the averaged 
entropy within the 10% of r^^^ ('S'o.ir2oo )> the ratio of the 
total internal to kinetic energies (E^/Ej^), and the inte- 
grated X-ray luminosity (Lx), both within the radius r^oo- 
As in previous plots, each column represents the results of 
a cluster representing one of the three classes. 

For the relaxed cluster, the situation is simple. The early 
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Figure 8. Fraction of all galaxy clusters in our simulation that 
have CC vs the gaseous mass of the clusters. The simula ted clus- 
ters (c ontinuous line) and the observed c l usters from Ch en et al.l 
(•2007*) (dashed line) and lO'Hara et al] ||2006|) (dotted line) are 
binned in five linearly equispaced bins. 



stages of cluster formation have left high entropy and in- 
ternal energy. However, with the time evolution, the clus- 
ter cools and loses internal energy, creating a lower entropy 
core, and increasing the X-ray luminosity. The minor merger 
cluster exhibits a different history. As the cluster forms a bit 
later than the previous one, it begins with a lower entropy 
and internal energy compared with the relaxed cluster. The 
time zone when mergers happen - delimited by the vertical 
lines - can be clearly connected with important changes in 
the cluster evolution. The first minor merger boosts the en- 
tropy level and the internal energy, indicating that some en- 
ergy has been injected in the system. This energy reheats the 
ICM and produces a decrease in the luminosity by delaying 
the cooling. Later on, the cooling takes over again dumping 
part of the energy, but leaving a net increase in the entropy. 
The history of the major merger cluster is slightly differ- 
ent. At the initial phase, the smooth accretion has produced 
an increasing trend in the core entropy, the internal energy 
and the luminosity. After the major merger the situation 
is different. Due to the more dramatic effects of the major 
mergers (higher disruption, stronger shock waves, and more 
turbulence and mixing) there is an increase in the entropy 
level, but associated with an immediate loss of energy due 
to radiation. Another minor merger produces some minor 
changes but the final state is a cluster pretty similar to the 
previous ones but with a significantly higher entropy. 



3.4 Cool cores and cluster mergers 

It is well known that clusters of galaxies exhibit an impor- 
tant feature that allows to classify them into two separate 



Galaxy cluster mergers 11 



CL03 CL10 CL12 




0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 




0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 




0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 



z z z 
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{Ejj) and kinetic (S^) energies, and the X-ray luminosity (i^)' both evaluated within Tjqq. Left column shows the result for a cluster 
representing the major merger clusters, central column for a minor merger cluster and right column for a relaxed cluster. Vertical lines 
delimit the time interval when mergers occur. 



populations, those having cool cores (CC) and those others 
not having cool cores (NCC). 

Recently, IChen et alj (|2007l ) concluded that, roughly, 
half of the population observed in a sample with more than 
hundred clusters have CCs. The explanation for this di- 
chotomy is not clear and it remains a matter of debate. Sev- 
eral authors have studie d this prob l em by means of numer- 
ical simulations. Thus, iKav et aP (|2007h overestimate the 
number of CC clusters since almost all their clusters show 



the presence of CCs. However, iBurns et al.l (120081 ) claimed 
to be the first authors producing a simulation with CC and 
NCC clusters in the same numerical volume, although their 
abundance at z ~ of CC cluste rs, ~ 16%, s e ems t o be lower 
than the observed fraction by IChen et al. I (l2007l ). ~ 46%. 
Interestingly, the results presented by Kay et a!Lr(l2007 ) are 
based on SPH simulations, whereas those of Burns et al.l 
(|2008l ) are obtained using an Eulerian AMR code. It is likely 
that feedback processes could be directly involved in the sur- 
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Figure 9. Fraction of CCs as function of rcdshift for the simu- 
lated clusters in our complete sample. Error bars show VlV un- 
certainties due to the number counts. 



vival of CCs in clusters, but it is also possible that mergers 
could play an important role erasing the presence of CCs 
l|Poole et al.l bOOa iBurns et allbOOSl ). In this Section, we 
analyse our simulation paying special attention to the pres- 
ence of CCs, a nd their r elative abundances. 

Following iBurns e t al. (200^, we define a CC cluster 
as one with a ^ 20% reduction of its central temperature 
compared with the surrounding region. Using this definition, 
we have classified all the clusters in our complete sample - an 
extended sample containing all the clusters in the simulation 
- into two groups: CCs and NCCs. 

In Figure [S] we plot the fraction of CCs as a function 
of gaseous mass at z = 0. We have binned the clusters using 
five linearly equispaced bins in the range 
The continuous line shows our results. For the sake of com - 
paris on, we have used the samples of lO'Hara et al.l (|2006l ) 
and IChen et al.l l|2007D . and we have binned the clusters in 
these samples using the same bins than for our results. The 
dashed and dotted lines correspond to Chen's and O'hara's 
results, respectively. 

Our re s ults a re extremely similar to those of the 
iBurns et al.l (|2008l ) simulation, where they found a total 
fraction of ~ 16% CC clusters (see figure 8 in this refer- 
ence) and i n our case the number is also ~ 16%. As in 
iBurns et all (|2008 



our results differ in the absolute num- 
bers from the observational data, but more interestingly, we 
have confirmed the general trend of a decreasing number of 
CC clusters with cluster mass. 

Although, we have no clear explanation for the discrep- 
ancy between the absolute number of CC clusters in our sim - 
ulation and the observational data by I Chen et al.l (|2007h . 
two plausible explanations can be given in order to interpret 
these results. The first one has to do with the fact that no 
metal-dependent cooling has been c onsidered in the simu- 
lation . It is known (see for in stance iDe Grandi fc Molendil 
l|2002l ); IVikhlinin et all (|2005M that some clusters can show 
strong metallicity gradients, with metallicities rising to solar 
in the central regions. This limitation of the present simu- 



lation could produce some artificial reduction of the cool- 
ing, specially in systems where kT < 2 keV. Therefore, this 
shortcoming could mimic, effectively, some sort of uncon- 
trolled non-gravitational feedback. The second possibility is 
related with a resolution issue, as no resimulations of the 
clusters have been performed. Therefore, despite the use of 
an AMR code, there could exist some resolution limitations. 
This last possibility seems much less important, though. 

We have also looked at the time evolution of the fraction 
of CC clusters. In Figure [O] we plot the fraction of CCs in 
our sample as a function of the redshift from z ~ 2 until z — 
0. Again, our results ar e fully consistent with the simulations 
by iBurns et al.l l|2008l ) and show no important change in the 
fraction of CCs backwards in time, at least back to z ~ 1. 
Our results are in contradiction with observational evidences 
showing an i mportant variatio n in the fraction of CCs from 
z = 0.5 ( Vik hlinin eral]|2006h . 

Before 2; ~ 1, we find a dramatic reduction in the frac- 
tion of CCs with time. As it would be expected, the abun- 
dance of CCs would be directly correlated with the hierar- 
chical formation of the clusters. At the epoch of cluster for- 
mation, almost none of the clusters would have a CC. The 
formation of CCs would require the establishment of cool- 
ing flows which, eventually, and through a slow process will 
form the cool cores. However, once the clusters were fully 
formed, the major mergers would destroy the CCs, creat- 
ing a population of NCC clusters. It is clear that feedback 
processes would also play a crucial role in this mechanism, 
but in the present simulation, where no relevant feedback 
mechanism - apart of the gravitational - has been taken 
into account, the effect of mergers on the existence of CCs 
is more outstanding. As the mergers are more dramatic in 
the more massive systems, this would explain the anticor- 
relation of the fraction of CCs and the mass of the clusters 
(see Fig.lH). 

In order to deepen our knowledge of the connection be- 
tween merger activity and the presence of CCs in clusters, 
we have studied the dependence of the fraction of relaxed 
clusters with the cluster mass. If mergers play an impor- 
tant role in the existence of CCs, one would expect that the 
systems that have evolved quietly (no recent mergers) do 
have cool cores. The comparison between the fractions of 
CC clusters and relaxed clusters is not direct as the estab- 
lishment of cooling flows and the subsequent formation of 
CCs could require long time scales, specially for the smaller 
systems. In any case, and taking the result with caution due 
to all the uncertainties, we have computed the fraction of 
relaxed clusters which have central cores with cooling times 
shorter or comparable to the elapsed time from the clusters 
formation {zfor, see Sec. 3) until the actual time {z — 0). 
Therefore, these clusters would have had time to set up a 
CC. The result shows the same trend than that in Fig. |8l 
that is, the number of relaxed cluster (no mergers) decreases 
with the cluster mass. This comparison would show that the 
smaller systems tend to have a CC and a quite evolution (no 
merger events), whereas the larger systems suffer the most 
important merger events and are NCC systems. 

3.5 Scaling relations 

The scaling relations are crucial tools to study the galaxy 
clusters, as they connect observables like X-ray luminosities. 
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Figure 10. Scaling relations for our galaxy cluster sample. The 
panels display, from top to bottom: integrated X-ray luminosity, 
mass, and mean entropy within the radius r^gg . In all the panels, 
the previous quantities are plotted against the temperature, T^qq , 
computed according to Eq. [9] The different symbols represent the 
clusters in the sample, the continuous lines stand for the proper 
fitt ings, and the sma ll filled dots represent observational dat a 
by iHornej l l200ll) and [P onman. Sanderson fc FinoguenOTlll2003l) . 
respectively. 



with cluster properties, namely, masses and temperatures. 
Moreover, they can be an excellent way to check the be- 
haviour and consistency of the simulations, by comparing 
with the scaling relations obtained in other simulations or 
with observations. 

The galaxy cluster reduced sample studied in this pa- 
per is biased towards the most massive clusters of our sim- 
ulation. Therefore, the statistical properties of this sample 
must be taken with caution as the sample is far from being 
complete, due to the numerical limitations. In the present 
subsection, we have extended the reduced sample (see Table 
1) by considering all clusters, in the complete sample, with 
temperature, T^^g ^ 1 keV (see Eq. [9]). 

In Figure [TOl several scaling relations are plotted: X-ray 
luminosity (upper panel), mass (middle panel), and mean 
entropy (bottom panel) within the radius r^„g. All these 
three quantities are plotted against the temperature T^og. 
Our results can be fitted by the following scaling relations: 
L cc M oc and S oc T°-^. In the three rela- 

tions, we have plotted all the clusters with temperatures 
Tjoo £ [1.0, 8.0] keV. This choice slightly increases the num- 
ber of clusters of the original sample presented in Table 1. 
For completeness, and in order to compare with observa- 
tional dat a, we h a ve com pared our scaling relations with 
data by iHorneil ll200l|) for the L — T relation and by 
IPonman. Sanderson fc Finoguenovl l|200^ ) for the S — T re- 
lation, respectively. These data are displayed as small dots 
in the top and bottom panels in Fig. 1101 The results for 
our simulated sample seem to be consistent with observa- 
tional data, leaving aside all the uncertainties of such direct 
comparison. 

Focusing on the efi'ect of mergers, and for the sake of 
comparison with previous works, let us assume that clusters 
that have had a relatively quiet evolution would likely de- 
velop a CC, whereas those clusters involved in merger events 
would see their cool cores distorted, turning into NCC clus- 
ters (see Sec. 3.4). Under this assumption, we could consider 
that the galaxy clusters in our sample, labelled as major and 
minor mergers, could be identified, bro adly speaking, wit h 
NCC clusters like the ones studied by iPoole et al] (|2008t ). 
and analogously, our relaxed clusters would be the CC clus- 
ters in this reference. The results displayed in Fig. [TO] show 
some degree of segregation, with most of the major and mi- 
nor merger clusters located at well separated regions on the 
scaling relation plots. The minor merger clusters sit, prefer- 
entially, at an intermediate region between the major merger 
clusters and the relaxed clusters. As a gross trend, the ma- 
jority of clusters which have suffered mergers, are placed in 
zones with higher temperature and higher luminosity, mass, 
and entropy, respectively. 

Mergers, specially the major ones, typically boost clus- 
ters along the L-T relation but not parallel to this relation. 
In Figure [TT] we plot - for the major merger clusters in Ta- 
ble 1 - the overall drifts experienced by these clusters from 
z — 1.5 until z = 0. This evolution is illustrated by vectors 
starting (finishing) at the values for L and T ai z — 1.5 
{z = 0). All the vectors could be decomposed in two com- 
ponents representing the change in temperature and lumi- 
nosity. Two special clusters deserve a particular discussion. 
The first one, represented by a triangle, has no arrow as- 
sociated. This is due to the fact that this cluster has been 
classified as major merger but it marginally satisfies the 1:3 
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Figure 11. The X-ray luminosity vs the temperature, T^qq, 
within radius, r^pp. The plot shows, as vectors, the overall drifts 
experienced by the most massive major merger clusters in their 
positions at the L-T relation from z = 1.5 until 2 = 0. The clus- 
ters, at 2 = 0, are labelled by squares, except two particular cases, 
identified by a triangle and a star, discussed in the text. Dashed 
line stands for the fitting of the L-T relation for the whole sample. 



condition in the mass ratio. Therefore, it is a transition case 
between major and minor clusters, according to our defini- 
tion of merger. However, as we will discuss later, consistently 
with its evolution in the L-T plane, it behaves like the minor 
merger clusters. The second particular case, represented in 
the plot by a star, shows the evolution of a cluster with a 
extremely strong cooling flow at z — 1.5. A careful study 
of the radial temperature of this cluster at that time, shows 
an extremely relevant CC. During the evolution, the merger 
event substantially reduces the cooling flow, and disturbs the 
CC, although the cluster remains radiating. Interestingly, all 
clusters, except the one represented by the triangle, show a 
net increase in temperature. Concerning luminosity, letting 
aside the two particular cases just mentioned, the rest of 
clusters shows an increase in luminosity. We have performed 
a similar analysis for the minor merger clusters, and they do 
not show any clear trend and seem to behave very similar 
to the relaxed clusters. 

There is a clear bias in the treatment of small clusters 
and groups. Nevertheless the results of our simulation are 
consistent with previous results, specially considering that 
it has not been introduced any other pre-heating or feedback 
mechanism, besides the one from the star formation (very 
poor in the present simulation) and the so-called gravita- 
tional feedback (shock waves, mergers, etc). In any case, 
even though the sample can be limited, we wish to stress 
that the individual properties of each of the most massive 
clusters are well defined. 



4 DISCUSSION AND CONCLUSIONS 

We have presented the results of a hydro and dark matter 
simulation of a moderate size volume of the Universe in the 
framework of a concordance cosmological model. The simu- 
lation follows the evolution of gaseous and dark matter com- 
ponents. Other relevant processes, like heating and cooling 
for a primordial gas, have been also taken into account. 

The main idea of the present paper has been to study 
the role of galaxy cluster mergers as source of feedback of 
the ICM. The general picture of our simulation shows how 
mergers stirred the ICM by producing shocks and sound 
waves in scales comparable to the dimension of the cluster. 
These phenomena produce turbulence in the form of large 
eddies characterised by Reynolds numbers. Re ^ 10^ , which 
are values accessible to present-day numerical simulations 
dRicker fc Sarazinll200ll). 

Recent work by lAgertz et al.l ('2007') has shown, in some 
idealised situations, how the different rmmerical techniques, 
namely SPH and AMR, produce different results when de- 
scribing the formation and evolution of eddies (Kelvin- 
Helmholtz) instabilities. According to this last reference, 
there are no substantial differences in results obtained with 
both techniques if the evolution time is smaller than the 
characteristic time of formation for Kelvin-Helmholtz insta- 
bilities, tkh- This characteristic time can be estimated for 
the case of a cluster merger assuming similar densities for 
the clusters, a characteristic length scale of a few hundred 
of kpc, and a relative velocity of the order of one thousand 
km/s; with all this conditions, tkh turns to be of the order 
of a few hundreds of millions of years. 

In the particular case of galaxy cluster mergers, the typ- 
ical time scales of evolution after a merger are much larger 
than this tkh, and therefore, it is likely that the particular 
numerical scheme used to simulate such scenario could play 
some role in the results. This has motivated us to study this 
problem by means of an AMR technique, and as a comple- 
mentary work to the studies already published. Therefore, 
we have used an AMR Eulerian code specially designed for 
cosmological applications, with excellent capabilities to deal 
with strong gradients, shock waves, and low density regions. 
These kind of codes also describe properly some instabilities 
phenomena such as Kelvin-Helmholtz, Ray leigh- Taylor and, 
in general, turbulence and mixing processes (see the recent 
work by MitchcU ct al. (2009)). 

In order to study the effect of mergers as feedback 
source, we have extracted and followed the evolution of the 
galaxy cluster like halos in our simulation. These halos have 
been studied directly from the simulation and without any 
resimulating scheme. This has implied a limitation due to 
numerical restrictions. Since the numerical scheme has the 
tendency to better resolve the most massive halos, our sam- 
ple is biased towards these systems. However, it is in those 
large systems where the effects of mergers would be more 
relevant. On the other hand, this apparently drawback of 
lack of resolution on small halos, is compensated by the ad- 
vantage that mergers can be followed in a consistent way as 
they naturally occur in the evolution of the simulated vol- 
ume of the Universe. Thus, no controlled collisions have to 
be imposed a priori. 

We have assumed a definition of galaxy cluster merger 
depending on the mass ratio of the halos involved in the 
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merger. As additional condition, we have introduced a time 
limitation, in such a way that only mergers occurring in the 
recent past are taken into account. Thus, if masses are sim- 
ilar - between 1 and 1/3 - we define those events as major 
mergers. Events with larger ratios in masses are classified 
as minor mergers. Finally, clusters evolving without rele- 
vant merger events are designed as relaxed. Once we have 
assumed this criterion to group the clusters in our simula- 
tion, we have studied the radial profiles of each cluster for 
the gas and dark matter densities, temperature and entropy. 
In order to compare the main differences among the three 
classes, we have computed average profiles for each group. 
The forms of the different profiles are basically the same for 
the three categories, indicating no substantial changes in the 
physics of clusters. However, there is a trend in the normal- 
isation. The relaxed and minor merger clusters have similar 
values of all quantities, whereas the major merger clusters 
are slightly hotter and with higher entropy. 

The previously mentioned trends in entropy and tem- 
perature of clusters depending on their evolutionary history, 
can be quantify by looking at a representative cluster of each 
class and comparing the time evolution of their profiles of 
temperature and entropy. 

In the same manner, the effects associated with cluster 
mergers can be traced in the time evolution of global quanti- 
ties like the entropy in the inner 10% of the virial radius, the 
ratio of internal to kinetic energy, or the X-ray luminosity. 

In all cases, cluster mergers release energy which ends 
up partially in the final object. The amount of energy locked 
in the final cluster is significantly larger for cluster mergers 
of similar masses (major mergers). 

We have also considered all the clusters in the 
simulation without differentiating amid the merging ac- 
tivities. These res ults have been c ompared with previ- 
ous simulations (IVoit. Kav fc Bryan 2005; iBurns et al.l 



I2OO8I: Kav et al. _ 20071) and with ob se rvational data 
llPonman. Sanderson fc Finoeueno^ l2003l : IVikhlinin et al.l 
I2OO5I : IPratt et al.ll2007l l. paying special attention to the en- 



tropy and temperature profiles. Our results seem to be con- 
sistent with both, simulations and observations, in an aver- 
age sense. However, there are still important differences like 
the lack of central drops in the temperature profiles of the 
most massive clusters. 

The fraction of clusters in our sample that has cool cores 
has been computed at several redshifts. At z = 0, our re- 
sults are fully comp atible with previous AMR simulations by 
iBurns et all (|2008l l. although seem to differ with the results 



of the SPH simulation by iKav et all (j2007l ). We have com- 
pared the fraction of c ool cores in o ur sim ulation with the 
observational data by IChen et al.l (|2007h showing a simi- 
lar trend, that is, the number of clusters with cool cores 
decreases with the cluster mass. 

Given the fact that in our simulation the gravitational 
feedback is the relevant feedback mechanism, we have tried 
to correlate the merger events with the existence of cool 
cores. In order to do so, we have computed the fraction of 
relaxed clusters (no merger activity) as a function of the 
cluster mass. Interestingly, the fraction of clusters with cool 
cores, and the fraction of relaxed clusters as function of the 
cluster mass, show a very similar trend. Unfortunately, the 
absolute numbers of CCs in our simulation and the obser- 
vations are quite different. We suggest two possibilities ex- 



plaining this discrepancy. The first one would be related with 
the fact that no metal-dependent cooling has been consid- 
ered in the simulation. This simplification could make the 
cooling more inefficient, specially at the central regions of 
the clusters. The second reason would be linked with a lack 
of resolution, which appears to be quite unlikely given the 
actual features of the considered simulation. In any case, it 
seems clear that there is an evident link between the merger 
events and the no existence of cool cores. 

On the other hand, the time evolution of the fraction 
of cool cores shows that this quantity has not changed sub- 
stantially from jz ~ to z ^ 1. This result i s compatible 
with previous simulations l|Burns et al.ll2008|) but in di s- 
agreement with observational data I Vikhlinin et al.ll2006l ). 



The cluster sample analysed in this paper is limited 
due to the fact that no resimulations have been done, and 
therefore, although an AMR code has been used, there are 
still resolution limitations. Nevertheless, we have analysed 
the scaling relations derived from our sample. Our results 
for L oc T^-^, M oc ^ and 5" cx T"'^ are consistent with 
previous results that do not introduce any extra reheating 
or feedback. 

We have found some degree of segregation in the scal- 
ing relations depending on whether the clusters have or have 
not undergone a recent merger. The systems that have ex- 
perienced merger events are usually located at high temper- 
atures, luminosities, masses, and entropies, respectively, at 
the different scaling relation plots. These results could be 
comparable with recent wor ks looking at the existence or 
not o f cool cores in clusters (|Burns et al ] |2008l : |Poole" et al.l 
l2008h . The analysis of the time evolution of major merger 
clusters in the L-T relation, has shown that these clusters 
have a tendency to move towards regions of this relation 
with higher temperature and luminosity . This tendency is 
similar to that found by iHartlev et al.l l|2008l ). where au- 
thors investigate the L-T relation in a large simulation with 
a strong preheating. 

A clear improvement for future work would be to in- 
crease the number of clusters in the sample by simulating 
larger volumes with higher resolution. Therefore, it would 
be feasible to reliably study the scaling relations for each 
one of the three families of clusters that we have considered 
in this paper. In any case, even when the sample can be lim- 
ited, the individual properties of each of the most massive 
clusters are well defined. 

The role of mergers as source of feedback, transferring 
part of the gravitational energy to the thermal energy, is 
still a matter of debate and study. Mergers are crucial to 
understand galaxy cluster formation and galaxy formation 
scenarios as they infiuence directly the ICM properties. Sim- 
ulations with higher resolution and including more physical 
processes are needed in order to keep on quantifying the 
role of mergers in the hierarchical scenario of structure for- 
mation. In parallel, some results from simulations, like the 
ones presented here, can be considered in semi-analytical 
models in order to improve their description of the gas com- 
ponent. 
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APPENDIX A: STAR FORMATION 

The star formation has been in troduced in the M AS- 
CLET code following the ideas of lYepes et al.1 (|l997l l and 
ISprineel fc Hernauisd (|2003l ). In our particular implementa- 
tion, we assume that cold gas in a cell is transformed into 
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star particles on a characteristic time scale t, according to 
the following expression: 

dp* dp 
~dt ^ 



dt 



= (1-/?) 



i.(p) 



(11) 



where p and p, are the gas and star densities, respectively. 
The parameter (3 stands for the mass fraction of massive 
stars (> 8Mq) that explode as supernovae, and therefore 
return to the gas component in the cells. We have adopted 
/3 = 0.1, a value compatible with a Salpeter IMF. For the 
characteristic star formation time, we make the common as- 
sumption = iS(p/Pth)~^^^, equivalent to p, = p^'^ /tl 

(jKennicutd l 1998:). In this way, we have introduced a depen- 
dence on the local dynamical time of the gas and two param- 
eters , the density threshold for star formation (Pj^) and the 
corresponding characteristic time scale {tg). In our simula- 
tions we have taken to = 2 Gyr and p^^ = 2 x 10~^^ g cm~^ . 
From the energetic point of view, we consider that each su- 
pernova dumps in the original cell 10^^ erg of thermal en- 
ergy. 

In the practical implementation, we have assumed that 
star formation occurs once every global time step, Afi=o 
and, only in the cells at the highest level of refinement {I = 
7). Those cells at this level of refinement, where the gas 
temperature drops below T < 2 x 10* K, and the gas density 
is p > Pj^ = 2 X 10~^^ (? cm"'', are suitable to form stars. 
In these cells, coUisionless star particles with mass m* = 
p»/S.ti^(,/S.x^ are formed. In order to avoid sudden changes 
in the gas density, an extra condition restricts the mass of 



the star particles to be m, — min(m, 
is the total gas mass in the considered cell. 



, 2 ' it'gas 



), where rag 



